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ABSTRACT 

We present the first astronomical observations obtained with an Apodizing 
Phase Plate (APP). The plate is designed to suppress the stellar diffraction pat- 
tern by 5 magnitudes from 2 — 9X/D over a 180° region. Stellar images were 
obtained in the M' band (Ac = 4.85/im) at the MMTO 6.5m telescope, with 
adaptive wavefront correction made with a deformable secondary mirror designed 
for low thermal background observations. The measured PSF shows a halo in- 
tensity of 0.1% of the stellar peak at 2\/D (0.36 arcsec), tapering off as r~^/^ 
out to radius 9\/D. Such a profile is consistent with residual errors predicted for 
servo lag in the AO system. 

We project a 5a contrast limit, set by residual atmospheric fluctuations, of 
10.2 magnitudes at 0.36 arcsec separation for a one hour exposure. This can be 
realised if static and quasi-static aberrations are removed by differential imaging, 
and is close to the sensitivity level set by thermal background photon noise for 
target stars with M' > 3. The advantage of using the phase plate is the removal of 
speckle noise caused by the residuals in the diffraction pattern that remain after 
PSF subtraction. The APP gives higher sensitivity over the range 2 — 5X/D 
compared to direct imaging techniques. 

Subject headings: Instrumentation: high angular resolution — stars: individual 
{/J, Her A) — stars: low mass, brown dwarfs 



^Observations reported here were obtained at the MMT Observatory, a joint facihty of the University of 
Arizona and the Smithsonian Institution. 
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Introduction 



The direct imaging of extrasolar planets (ESPs) presents a daunting technical challenge. 
The small angular separation from the host star places the planet within a diffracted starlight 
halo many decades brighter than the planet itself. A coronagraph suppresses this diffraction 
whilst preserving the angular resolution and flux of the planet through to the final imaging 
camera. There are many coronagraphic designs and they all represent the fundamental 
trade-off between throughput, angular resolution and diffraction suppression as a function 
of position in the final image. 

The prototypical coronagraph jLvotlll939h is an imaging camera consisting of optics that 
form two intermediate focal planes. The first is an image plane where an occulting mask is 
used to remove the Airy disk and diffraction rings of the central star. A sharp edged mask 
in the image plane introduces high spatial frequency components in the wavefront, which 
are masked out in the subsequent pupil plane with an undersized pupil stop (the Lyot stop). 
The original design is not optimal for ESP detection, and results in reduced throughput and 
angular resolution. 

Encouraged by t he discovery of near ly two hundred ESPs through radial velocity searches 
around nearby stars (IButler et al.ll2006l ). the potential of directly imaging terrestrial sized 
planets has fueled research in a wide range of coronagraphic designs, a large selection of 



which are reviewed, analyzed and compared in iGuyon et al.l (120061 ). In that paper, they 



concentrate on designs suitable for space-based telescopes consisting of an unobstructed cir- 
cular aperture, which produce contrast ratios in the regime of 10~^°. Most of the designs 
considered modify the incoming wavefront by placing apodizing masks in the image plane 
and/or in the pupil plane of the coronagraph. The designs can be further subdivided by 
whether they modify the amplitude or the phase of the wavefront, or some combination of 
the two. 

Most of the image plane apodizing me thods are sensitive to tip-tilt errors and the finite 
size of the stellar disk (IGuyon et al.ll2006l ). allowing light from the primary star to add an 
extra source of noise and reduce the achievable contrast (iKuchner fc Traubl |2002| ). Recent 
theoretical studies have indicated that eighth-order band limit ed masks, which are designed 



to be insensitive to srn all amounts of low order aberrations (IKuchner. Crepp &: Gd 12005 



Shaklan fc GreenI |2005| ) , show promise for providing deep suppression over small opening 
angles and with theoretical throughputs of up to 50% and IWAs (inner working angles) of 
4A/-D. Laboratory experiments on manufacturing proce ss limited masks ga ve a measurement 
of 15% throughput for a designed transmission of 20% (ICrepp et al.ll2006l ). 



Alternative methods for suppressing diffraction rely on variable transmission masks 
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in the pupil plane, either as binary mas ks (IKasdin et al 



2003 



as graded (apodized) transmission masks (INisenson fc Papaliolios 



Vanderbei et al.l 120031) or 



2001 



Aimell2005f ). Pupil 



apodizing methods generate field invariant PSFs and so do not suffer from the tight alignment 
tolerances required by image plane masks, and their suppression is not affected by the angular 
size of the target star. 



Guyon et al.l (120061 ) considered a phase apodization algorithm from I Yang fc Kostinski 



(I2OO4J ). which uses a square pupil and provides deep suppression over a 90 degree opening 
angle with a 60% loss of light from the central PSF core. Our independently developed phase 
apodization algorithm can generate PSFs for more generalized pupil geometries (including 
pupils with secondary obscurations) over an opening angle of 180 degrees. For the phase 
plate in this paper, the FWHM of the Airy core increases by 10% and produces a considerably 
smoothed zone of suppressed diffraction from 5 to 6 decades below that of the core, at a 
modest cost of light (29%) from the central Airy core of the imaged object and an IWA of 
2.2\/D. 

The sensitivity of close-in searches is increased by the subtraction of the PSF of the 
telescope, determined either by an analytical model, images of a nearby reference star of 
similar brightness, or by using the rotation of the target star on the night sky to construct a 
model PSF with the putative companion signal removed by median combining the individual 
exposures. In the ideal high signal-to-noise limit, the PSF images subtract cleanly, leaving an 
image limited by photon noise following the diffraction pattern. In reality, PSF subtracted 
images show considerable small scale structure, and these systematic aberrations (thought 
to be caused by optical flexure within the telescope and camera system) limit the detection 
sensitivity of faint companions at small angular separations. 



The 6.5m MMT AO system (IWildi et al.ll2003l : iBrusa et al.lll999r i is the flrst reahzation 
of a large aperture telescope coupled with a deformable secondary mirror that delivers an AO- 
corrected //15 beam. At longer wavelengths (> 2.5/im), the MMT AO system is uniquely 
sensitive because of the lower level of background light it emits at infrared wavelengths 
compared to conventional AO systems on larger telescopes (iLloyd-HartI |2000| ) - typically 
twice as efficient as the Keck telescopes at 5yum. Furthermore, the mid-infrared wavelengths 
are well suited to the detection of hot extrasolar planets due to their enhanced thermal 
emission detectable through the E a rth's atmospheric transmissiori win dow at 5 microns 
( iBurrows. Sudarsky fc Lunind l2003l : iBurrows. Sudarsky fc Hubenyl 120041 ) . Where thermal 
telescope emission is the dominant source of background photons, a reduction in emissivity 
corresponds to a direct reduction in exposure time to reach a particular sensitivity. Using 
measurements of sky frame fiuxes at two differing airm asses, the emissivity of th e telescope, 
AO system and Clio camera is estimated to be 10% (ISivanandam et al.l l2006l ) . This is a 
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significant improvement over the 25-50% of other conventional mid-IR AO systems. 

In this paper we present a new method to suppress diffraction in an astronomical image: 
the Apodizing Phase Plate (APP). In Section [2] we present the basic theory behind the optic, 
and Section [3] describes the first on-sky APP observations and the data reduction. Section H] 
compares the expected point spread function using the APP with the MMT AO system, and 
the calculated sensitivity and limits to the achieved contrast ratios are in Section [51 Section 
E] presents an initial scientific observation using the APP and we discuss the implications for 
future observations in Section [71 



2. Phase Plate Design and Manufacture 

The starlight halo against which detection of faint companions are to be made consists 
of a deterministic diffraction pattern (or point spread function, PSF), and a speckled halo 
caused mostly by residual phase errors left over by the AO system, and to a lesser degree 
by scintillation over the pupil. There are also inevitably going to be optical imperfections in 
the optics, e.g. small phase aberrations and transmission variations that will lead to faint 
static and semi-static speckles in the halo. For the MMT, these have not yet been seen or 
characterized. The scintillation halo is expected to be fainter than the AO-residual phase 
halo, and at high Strehl, the diffraction pattern is the brightest. Reducing the diffraction 
halo is the first step in improving the S/N of a faint planet near the star, since it reduces 
the "noise" against which the detection is made. Preserving the light in the diffraction core 
is also important since it affects the "signal" coming from the planet. The phase plate is 
intended to suppress the deterministic diffraction pattern to a level below that of the next 
limiting component of the halo (i.e. the residual AO phase errors), while doing the least 
harm to the light in the diffraction core. 



he phase plate was designed using the approach outlined as "Method I" in lCodona fc Angel 



( 2004). This method has recent ly been validated in a closed-loop laboratory experiment 



([Putnam. Codona fc Angell l2006l ) . The underlying principle is to create an "anti-speckle" 
for each "speckle" or diffraction structure within some region of interest (ROI) within the 
PSF of the star. (An "anti-speckle" is an intentionally-constructed region of the halo which 
matches an unwanted speckles' location and amplitude, but is 180 degrees out of phase. The 
anti-speckle is intended to cancel a given speckle). We accomplish this by introducing a weak 
sinusoidal phase "ripple" across the pupil which acts as a weak diffraction grating, causing a 
small amount of the starlight core to be diffracted into a new halo speckle. Strictly speaking, 
a line of speckles are formed corresponding to the higher diffraction orders of a grating, but 
in the limit of a small sinusoidal ripples, the higher orders can be neglected. By using a 



- 5 - 



phase ripple with spatial frequency k = 27r6/X, we can create a speckle at an angle 6 relative 
to the star. For a small phase ripple, the amplitude of the created speckle is proportional 
to the ripple amplitude, and the speckle phase relative to the core is linearly related to the 
ripple phase. Since the phase ripples are real (i.e. only affect phase and not amplitude), each 
created speckle has an anti-Hermitian (antisymmetric real, symmetric imaginary) counter- 
part on the other side of the star. This is the opposite symmetry of diffraction patterns, 
and therefore it is only possible to suppress diffraction over at most one half of the region 
around the star. The region on the other side of the star is reinforced, adding energy to the 
diffraction pattern there. 

By linearly superimposing a set of such ripples, we are able to reduce the halo over 
the region of interest (ROI). Since this algorithm does not take into account the detailed 
shape of the created speckles, and the process of creating them alters the initial halo in a 
nonlinear fashion, this approach will only suppress, not cancel, the halo in a single step. 
However, by iteration, we can robustly find our way to a solution which nearly cancels the 
ha lo over the ROI- For small iteration steps, the formula for the n + 1 iteration of the phase 



is f jCodona et all 120061 ) 



where Il{x) describes the pupil, Ti.{x) = ^^qi^^'^ ^-^P (^'^ ' ^) /(^^)^ spatial filter 

corresponding to the dark region of interest in the focal plane, and takes the imaginary 
part of the result. Equation [T] is derived by making the halo within the ROI dark without 
regard for the effect on the PSF core or the halo outside of the ROI. While this approach 
worked adequately in this case, it is sub-optimal in terms of Strehl ratio. The theory has 
now been extended to include the maximum preservation of power in the diffraction core. 
The extended theory will be presented in a later paper. 

Suppressing diffraction by introducing aberrations is inevitably chromatic. The prac- 
tical effect is that the dark zone in the halo begins to brighten with increasing bandwidth. 
However, so long as the combined halo remains at least an order of magnitude below the 
residual AO speckle halo, it won't significantly limit detection. Since the mismatch between 
the complex antihalo and the halo grows linearly for small variations in wavelength, the 
intensity of the halo floor brightens like (A — A n)^ where Aq is the wavelength where the eq. 



[T] iteration was performed (ICodona et al.ll2006l ). Integrating the growing halo over the entire 



detector band shows that the total halo background grows as the cube of the bandwidth. 
This leaves us a reasonable bandwidth with which to work, effectively controlled by the 
Strehl ratio achieved by the AO system. In general, depending on the science goal, working 
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with a wider bandwidth is better, since it collects more photons and gives better speckle 
noise averaging. Therefore, it is best to use as wide of a bandpass as possible until the 
diffraction halo starts to become significant. 

The mean residual AO speckle halo is unaffected by the phase plate and remains at the 
same level, while the Strehl ratio of the phase plate drops the PSFs of the star and planet 
into the noise. For this reason, it is better not to let the iteration run to something like 
convergence, but rather to stop it as soon as the solution becomes "good enough." Since the 
Strehl ratio starts at unity and drops monotonically as the iteration is performed, it is easy 
to place a threshold and stop the calculation when the desired level is reached. 

Two APPs were manufactured to our computed specifications by //- VI Incorporated, 
both using zinc selenide (ZnSe, n(4.68/im) = 2.43) as the transmitting substrate for the 
phase pattern. The wavefront surface was generated using a diamond tool whose pitch 
depth was synchronized with the rotation angle of the machine lathe holding the APP. The 
first APP used a phosphor bronze plate aligned with the phase plate pattern to mask the 
telescope pupil and secondary mirror support structures. Initial tests with this first APP 
validated the APP method, but the plate has no anti-reflection coatings on either side, 
decreasing the transmission efficiency of the optic and adding a scattered halo of light from 
the two internal reflections of the optic. The success of this plate demonstrated that pupils 
of arbitrary geometry, such as segmented mirrors on space telescopes, can be accommodated 
into the APP algorithm and produce a nulled region near the target of interest. We carried 
out scientific observations (detailed in Section E]) with this engineering version, and we were 
encouraged with the results to make a second, science grade APP. 

In the case of thermal imaging at the MMT, we also concluded that the thermal emission 
from the secondary support vanes does not significantly contribute to the noise background 
of our observations. Removing the secondary support masks from the APP design has two 
advantages - (i) the APP does not have to rotate to keep aligned with the secondary support 
vanes, resulting in cleaner background subtraction at thermal wavelengths, and (ii) the APP 
surface does not contain high spatial frequency regions required to deal with the presence of 
the secondary support vanes, resulting in a wavefront that is less sensitive to manufacturing 
errors in the diamond turning process. 

The science grade APP is manufactured as a single optic that contains the phase plate 
pattern, with a gold layer coating that blocks light outside the telescope pupil and within the 
secondary mirror's shadow. A broadband coating of magnesium fiouride with a measured 
reflectivity of less than 0.5% (from 3.2 to 5.5/im) provides high optical throughput and 
minimizes ghost image reflections. Figure [1] shows the design of the APP and the resultant 
PSF. The plate is designed to suppress diffraction from 2 — 9X/D in the focal plane of the 
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camera over a 180° region, forming a D-shaped region of suppression next to the target star. 
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Fig. 1. — Phase map of the APP design and its resultant theoretical point spread function. 
The left hand image shows the phase plate design cut into the zinc selenide plate, with the 
scale bar on the right showing height in microns. The secondary obscuration is deliberately 
over sized to allow for rapid in-situ alignment. Shown on the right hand panel is the calcu- 
lated PSF logarithmically scaled over 6 decades normalized to the peak intensity. The pixel 
scale on the right hand panel is for 5 micron imaging at the MMTO 6.5m telescope. 



The surface accuracy of the APP was measured using a Zygo optical interferometer 
at the manufacturers optical laboratory. The mean departure from the model surface is 
measured to be 70 nm r.m.s. Using the Marcchal approximation for estimating Strehl ratio 
at the operating wavelength of 5 microns, this corresponds to approximately 0.7% of the 
incident flux scattered out of the transmitted beam due to microroughness in the APP. 

The inner radius is set by signal to noise considerations for faint companions. Suppres- 
sion of the Airy pattern comes at a cost of energy from the core, and for this design we chose 
an APP core signal (the flux enclosed within the first dark Airy ring) that was 69.4% of the 
direct imaging PSF core signal. The outer working radius of the nulled region was set to 
match the control radius of the MMT AO system. Beyond the first Airy ring, the suppressed 
halo is designed to be 5 {3X/D) to 6 {6X/D) magnitudes fainter than the direct imaging PSF 
at a similiar radius, although we did not actually validate that claim by measurement. The 
phase plate theory is capable of creating much deeper suppressions, although in practice the 
limit will be set by the microroughness of the manufacturing process. One notable advantage 
of using the APP at the pupil plane of the telescope is that the suppression is independent 
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of tip-tilt and pointing errors in the telescope, which allows easy beam switching on the 
infra-red detector. 

Although it is possible to use the deformable mirror of an AO system to produce the 
APP wavefront, it is more practical to implement a single optic that sits downstream of the 
AO control loop. The two main reasons for this are: (i) the prescribed wavefront has many 
high spatial frequency features that require a high degree of actuator control to hold in place, 
and (ii) the wavefront sensor camera will see the wavefront superimposed on the atmospheric 
turbulence, potentially saturating the sensor. Keeping the APP out of the control loop allows 
the deformable mirror to be used in the removal of low-level residual static aberrations in 
the ROI of the A PP, and to test any cand idate companions for coherence with the light from 



the primary star (IKenworthy et al.ll2006l ) 



Observations 



Clio is an imaging camera (iFreed et al.l |2004J : ISivanandam et al.l |2006| ) designed for 
obtaining high spatial resolution images with optimum efficiency at L' and M band. It 
is optimized for imaging extrasolar planets at wavelengths where they are expected to 



have f luxes significantly in excess of that of a similar temperature blackbody ( IBaraffe et al. 
20031 : iBurrows. Sudarsky fc Lunind l2003l ). The plate scale of Clio is determined to be 



0.048574 ± 0.000090 arcsec pixel"^ from observations of two binary systems (HD 100831 
and HD 115404) in April 2006. The plate scale error is dominated by uncertainty in the true 
separation of these binaries. Position angle (PA) calibration was also determined from the 
two double stars mentioned above, and the accuracy is 0.20 degrees, again hmited by the 
uncertainty in the true position angles of the double stars. 

The APP was placed in the Clio camera at the intermediate pupil plane and mechanically 
positioned into the telescope beam using a motorized filter wheel. The camera has a pupil 
imaging mode which allows quick and easy alignment of the APP with the telescope pupil. 
To go from direct imaging to APP imaging takes less than five minutes, which is in strong 
contrast with systems involving an occulting mask in an image plane that require precise 
alignment for each new target observation. 

The APP was used on M' (A^ = 4.85/im) observations of HD 213179 (V=5.778 SpT=K2II) 
taken on 2006 July 09 12:10 UT. Extrapolating to M' magnitudes using color tables in I Cox 
(I2OOOI ) this star has an approximate magnitude M' = 3.1, which is the typical brightness 
of stars looked at in extrasolar planet surveys carried out with the Clio camera. The data 
were taken at an airmass of 1.06 in conditions of variable seeing, typically 0.6 — 1.0 arcsec at 
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500nm as estimated by wavefront sensor camera residuals in the AO system. We obtained in 
a set of 16 5-second exposures on target and 16 off target for a total of 80 seconds on-source 
integration. 

Approximately 1% of the pixels on the detector show anomalous behavior, including high 
dark current, high or low sensitivity, and time- varying responsitivity. Out of the methods 
we explored, we found the most robust way to select anomalous pixels was by looking at 
the sky-illuminated pixels in the beamswitched images. For a set of exposures making up 
a single beam, the variance of the sky illuminated pixels is calculated. Any pixels that 
have significantly high or low variance are marked as bad pixels. Combining the bad pixel 
masks from the two beam pointings gives a bad pixel mask that is applied to all of the 
unshifted science frames from the camera, interpolating over the bad pixels using adjacent 
good values. A blank region of the array is then used to sample a residual offset caused by 
the time-varying sky background, which is then subtracted from the image. The cleaned 
images are co-aligned using the IRAF imalign routine, which performs a cross-correlation in 
a box of 5 X 5 pixels centered on the unsaturated target star to locate the centroid of the 
PSFs. Resampling of the images so that the centroids of the PSFs are coincident is done 
using a cubic polynomial interpolant. These are combined with rejection of the most extreme 
pixel values to form the final image. The cleaned and background subtracted images were 
coadded to yield the summed image in the upper left panel of Figure [2l The PSF with the 
APP can be compared with that for the direct telescope PSF in the same figure. The Airy 
rings are strongly suppressed in the APP image. 

The utility of the APP becomes apparent when compared to direct imaging and PSF 
subtracted images. Figure [2] shows the telescope PSF comparing direct imaging with APP 
imaging. The Airy rings are strongly suppressed in the APP image, with only the time 
smoothed servo lag halo visible within 0.5 arcseconds of the target star. This is very much 
in evidence in the PSF subtracted images in the lower row. Here, four fake sources have 
been added to half the dataset and the other half of the data are used to subtract off the 
PSF. The background noise has been scaled in both the direct and APP imaging to represent 
identical on-source integrations of 40 seconds through the M' filter. The APP is background 
noise limited up to 2\/ D of the target star, whilst in the direct imaging case the closest fake 
source is confusion noise limited by the Airy ring residuals. 

The "Measured Light Profile" line in the upper panel of Figure [3] represents the measured 
light profile of the Clio APP image for a total on-sky integration of 80 seconds, normalized 
to the mean flux in the Airy disk of the star. The flux at a given radius is calculated by 
azimuthally averaging over a 150 degree wedge centered on the target star in the middle of 
the ROI. This curve can be compared with the dotted line labeled "Turbulence-free Model" , 
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which represents the theoretical APP PSF in the hmit of perfect AO correction over the M' 
band. Comparing the measured hght profile and turbulence-free model in the upper panel 
of Figure [31 we see that the measured starlight in the region of the suppressed halo is 4 to 
5 magnitudes brighter than the theoretical profile. This brightening is due to the residual 
aberration of the AO corrected wavefront, as discussed in the next section. 



4. PSF modeled for the AO- corrected wavefront 

The Strehl ratio for our observations is typically around 90% (see next section for calcu- 
lation). The ultimate limit to high-contrast imaging is from low-order (a few cycles/aperture) 
wavefront errors which cause speckles within a few X/D of the star. The diffraction suppres- 
sion provided by the APP allows us for the first time to directly measure the residual speckle 
intensities resulting from low-order (a few cycles/aperture) wavefront errors, allowing us to 
estimate the low-order performance of the AO system. 

Using the APP has an effect on both the signal and the noise portions of the detection 
sensitivity. The APP is designed to improve detection sensitivity by suppressing the diffrac- 
tion halo, but at the cost of light taken by the plate from the PSF core. Furthermore, other 
sources of background light will become the sensitivity constraint once the diffraction halo 
has been sufficiently suppressed, making deep halo suppression and its concomitant loss of 
Strehl ratio a liability. For a bright star, the next constraint in the sensitivity is the residual 
AO speckle halo, while for a fainter star the constraint might be photon noise from the 
thermal sky background. The azimuthally averaged light profile for the measured Clio APP 
data is shown as the solid line in the upper panel of Figure El 

Over the suppressed-diffraction region of the halo, the residual scattered halo is con- 
sistent with an angula r dependenc e of 9~^^^, which can be explained with a time-lag in the 



applied AO correction (lAngelll2003l ). This halo results when the wind advects the aberration- 
causing turbulence past the telescope aperture and there is a processing lag, Tiag, between 
the wavefront measurement and the application of the corresponding correction. Since the 
processing lag and the wind velocity multiply to give the shift between the actual wavefront 
and the AO-corrected version, the lag error can be parameterized by the distance diag which 
is the length of VTiag where v is the mean wind velocity. 
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Fig. 2. — Comparison of AO imaging with and without the APP. The upper row shows 
AO+APP imaging of a single star (HD 213179) and AO imaging of a binary system [jiHerA), 
with the flux and sky background of both images scaled to represent identical exposure times 
of stars with identical magnitudes. Both figures are logarithmically scaled so that the APP 
peak intensity has 69% of the no plate imaging intensity, consistent with the APP design. 
The lower two panels show the detection sensitivity of the two PSFs under PSF subtraction. 
Fake sources with a magnitude 7.5 fainter than the star have been added at 0.4 arcsecond 
increments {2X/D) from 0.4 to 1.6 arcseconds. The APP image intensity is scaled so that 
the fake source has the same signal in both PSF subtracted images, showing that even with 
the lower signal to noise the sources are more easily identified with the smoother, suppressed 
background and are not lost in the systematic Airy ring PSF subtraction. 
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Fig. 3. — Measured light profile and contrast ratio of the APP. The upper panel shows the 
azimuthally averaged light profile for a 150 degree wedge centered on the PSF for different 
cases described in the text, and the lower curve shows the Icr noise level, representing the 
sensitivity curve for the upper plot. The procedure for generating these curves is discussed 
in Sections m & [5l 
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To understand and quantify the AO halo in conjunction with the APP, we simulat ed our 



AO sy stem with the APP phase profile to model the combined PSF. As described in lAngel 



( l2003l ). a single Taylor wind flow results in a halo of speckles that are not evenly distributed 
with azimuth angle but instead form a two-lobed distribution (proportional to cos^cj)) 
aligned with the wind direction projected on the sky. Since this was not an obvious feature 
of our data, we modeled the wind as two superimposed turbulence layers oriented 90° from 
each other, resulting in a reasonably isotropic speckle halo. The two si mulated wind s were 



given the same speeds, but with independent Kolmogorov phase screens ( lHardylll998l ). The 
instantaneous residual phase error was computed by shifting and summing the two wind layer 
phase screens to the positions appropriate to time t, and then subtracting a low-pass-filtered 
version of the combined phase appropriate to time t — Tiag. The resulting residual wavefront 
was then passed through a Fourier optics model of the MMT, including the APP phase, 
resulting in an instantaneous PSF of the telescope. The phase screens were computed using 
a Fourier synthesis method on a 2048 x 2048 point grid, modeling a periodic 32m square 
region characterized by Kolmogorov turbulence with a characteristic Fried length tq. The 
evolving PSF was computed for 256 steps across the screen, shifting the array to simulate 
the wind evolution over a time of 32 /v seconds. Each frame represents an instantaneous 
realization of the atmosphere. By averaging together all the 256 frames, the resultant PSF 
represents a telescope integration that covers a 32m patch of atmosphere moving across the 
telescope pupil. The 5 second exposures from the Clio camera fix the exposure time, and so 
a PSF for a given length of atmosphere VTexposure (and therefore a given wind velocity) can 
be constructed by adding up a smaller subset of the 256 frames. Since a smaller number of 
exposures sums over fewer speckles, the variance in the individual exposures will be higher. 

By varying the wind lag distance and the Fried length in the AO models, we produced 
different modeled contrast ratio curves, the best fit of which is shown as the "Turbulence 
Model" in the upper panel of Figure [31 The best fit parameters are tq = 20 cm at the 
wavelength of 500 nm and a wind lag distance of 19 cm in each of the two wind directions. 

These parameters were then used to generate 16 simulated exposures, processed in an 
identical manner to the real dataset. Both the 16 Clio exposures and the 16 simulated 
exposures are averaged and the per pixel variance of the individual exposures about the 
resultant mean was calculated. This variance is divided by the number of images averaged 
together and then square rooted. The azimuthal mean of this quantity is then plotted as a 
function of radius over the image region identical to that in the upper panel, and we refer to 
these as our "Measured Sensitivity Curve" and "Turbulence Model Sensitivity for 80 seconds 
exposure" in the lower panel of Figured These represent the la sensitivity limit of the APP 
in detecting a faint companion using model PSF subtraction (see next section). 
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The width of the sensitivity curve is determined by the wind speed during the 80 seconds 
of total integration. By varying the number of summed frames from the total block of 256 
frames, the best match to the sensitivity curve was determined to be 192 frames, implying 
that (192/256)(32m) ^ v x 0.5s, which yields a mean wind velocity oi v = 4.8ms~^. Recent 



measurements of wind velocity as a function of altitude using SCIDAR (jPrieur et al.ll2004l ) 
show that the velocity of turbulent layers at altitudes of 10km can be considerably larger than 
the wind velocity measured at ground level (25ms~^ compared to 3ms~^ for observations 
at San Pedro Martir in 2000), and this trend is consistent with our measured ground layer 
wind velocity of 1.4 — 2.5ms~^. 
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Fig. 4. — The enclosed flux as a function of radius for various PSFs. This plot compares the 
direct imaging PSF with the idealized APP PSF and the measured APP PSF. The radius 
of the first dark ring indicates the flux enclosed within the central Airy disk. 

To determine how well our AO system was performing with the APP at 5 microns, we 
wanted to determine what the Strehl ratio of our data were. Figure H] shows the the fraction 
of the total flux enclosed within a given radius for various PSFs. The energy within the 
radius of the first dark ring indicates the fraction of energy in the image core. The dotted 
line represents the PSF of the telescope pupil, showing that the encircled energy of the core 
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is 74%, consistent with a pupil with central obscuration. The dot-dashed curve represents 
the modeled PSF of the telescope with an ideal APP, which is 69.8% of the encircled energy 
of the ideal PSF, or 52%. The data from Clio with the APP gives an encircled energy of 
45%. The ratio of the encircled energy from the measured APP to the ideal APP gives an 
estimate of how efficiently the AO system is working at keeping flux within the core of the 
image (and hence, the flux of the planet in its core) and this is approximately 86%. 

This encircled energy ratio is analogous to the Strehl ratio used in quantifying the 
wavefront correction of AO systems, and this is consistent with the Strehl ratio expected for 
the MMT AO system running at 550Hz with a 56 Zernike mode level of correction. Using 
the Marechal approximation we derive a mean wavefront error of 220nm. 



5. APP Sensitivity Comparison 

An observation looking for point-like sources goes further by performing PSF subtraction 
on the data. For scientiflc observations over many hours of on-source exposure, the camera 
and phase plate are flxed relative to the telescope during the set of exposures, so that the sky 
and background light distribution within the whole optical system does not systematically 
change. For our alt-az mounted telescope, the sky then rotates as a function of hour angle. 

We construct a master PSF by taking the median of the individual exposures, which 
rejects the flux from faint compansions whilst building a high signal to noise PSF that 
includes the static aberrations in the system. This master PSF is subtracted off the individual 
exposures, and each frame is rotated in software to align the sky orientation with the image 
orientation. The images are added together to form a flnal PSF-subtracted image. 

The "Measured Sensitivity" and "Turbulence Model Sensitivity" in the lower panel of 
Figure [3] closely match each other, and so we conclude that the total integration of 80 seconds 
represents a significant sampling of the turbulence power spectrum of the atmosphere. If we 
assume that the sensitivity curve will increase as the square root of the number of coadded 
exposures, we can estimate the expected sensitivity for a one hour exposure, labeled as 
"Turbulence Model Sensitivity for 1 hour integration". This sensitivity curve is dominated 
by two compo nents - one attributed to the photon shot noise of the sky background which 



scales as t^'^ (IHinz et al.l l2006l ) and is present at large radn, and a residual speckle halo 
noise component. In the limit of the residual speckle halo noise and sky background noise 
statistics both reducing as the halo noise is fixed with respect to the sky background 
noise, and the transition point for background noise dominating over halo noise is at about 
0.8 arcseconds {AX/D). If the star is brighter than M'=3, then the residual speckle halo 



-16- 



noise will be dominant out to larger radii, and if the star is fainter, the sky background noise 
will become dominant at smaller and smaller radii. 

Assuming that other static and quasi-static aberrations do not become dominant in 
longer integrations, we expect to be able to detect sources of 5M' — 10.2 with a Scr confidence 
level in one hour at a separation of 2\/D — 0.36" with the estimates of sensitivity made 
with our AO simulations. Our next investigation will be to obtain a much longer exposures 
with the APP and determine what sources of noise appear as the sky background noise is 
reduced, and to develop a suitable theory to characterize them. 




-2 2 

Fig. 5. — Image of the ji Her A system taken with an earlier prototype APP. Note how the 
secondary companion (to the lower right) has the same PSF as the primary star. 
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6. The fi Her A System 



In order to demonstrate that the principle of the APP was vahd, we observed a recently 
nearby binary sy stem, with the purpose of confirm ing the nature of the fainter component 
as determined by lDebes. Ge &: Chakrabortyl (|2002| ) and to obtain astrometric measurements 
for determining the dynamical mass of the system. 

The star /i Her A was observed with the earlier prototype of the APP at 2006 April 13 
UT 12:48 using the Clio camera at the //15 Cassegrain focus of the MMTO 6.5m telescope 
on Mount Hopkins, with the AO system operating in closed loop on the target star. The 
data were observed in the M' filter (Ac = 4.67/im, Xfwhm = 0.25). 

Our astrometric and photometric observations of the fi Her A system are listed in Table 
[1] with the image shown in Figure [51 Although we performed relative photometry in both 
M and M', our images are saturated (and thus non-linear) by a few percent in the M band 
images, and so that relative photometry is not reported here. iKidger fc Martin-Luid (120041 ) 
give J, H, K photometry for the /i Her A system that used a 15 arcsecond aperture, including 
both the components. Our M' band differential photometry, combined with the typical 
{K — M = 0.005) color of fi Her A and mK{Aa) = 1.743 results in an absolute magnitude 
of = 7.13 ± 0.05. The resultant colors for /i Her Ab are {H - K) = +0.6 ± 0.12 
and {K — M') = +1.17 ± 0.12, where errors are predominantly from the H and K band 
magnitudes, Mk = 8.3 ± 0.1 and Mh = 8.9 ± 0.1. 



Ex amining the p h otoni etry of 3 Gyr M, L and T dwarfs presented in iGolimowski et al. 



( I2OO4J ): iLeggett et al.l (120021 ) . we find that the absolute magnitudes are in dicative of a mid- 



M dw arf, consistent with the [H — K) color. Using the polynomial fit in IGolimowski et al. 
( I2OO4J ) for the spectr al type as a function of Mm', [i Her Ab is determined to be M4±l. Using 



the photometry from lReid fc Cruzl (120021 ) the Mk implies a later spectral type, around M6-8. 
Even assuming a later type star, the {K — M') color is still a whole magnitude too red for 
an M4 dwarf but similar to that for an early T dwarf, a conclusion ruled out by the age and 
the absolute magnitude of the object. Low resolution infra-red spectroscopy of this low mass 
object, along with further astrometric observations over the next decade will help determine 
its true nature. 



7. Conclusions and Future Plans 

We have demonstrated for the first time that an apodizing phase plate, generating a 
highly non-symmetric wavefront, can be manufactured and used in an imaging camera and 
telescope with AO system and produce a scientifically useful degree of diffraction suppres- 



-18- 



sion. The APP shape can also be introduced using a deformable element with high spatial 
frequency in upcoming AO systems to null out regions of interest, but this requires a high 
degree of closed loop control and also limits the bandwidth available for atmospheric cor- 
rection. The APP wavefront has regions of locally large wavefront tilt, and if introduced 
into an AO loop, it can cause Shack-Hartmann wavefront sensors to become saturated. By 
using a simple optic downstream of the wavefront sensing and correction of the AO system, 
we eliminate the need for these schemes and can implement, test and refine APP designs in 
short order. 

The APP gives a significant detection advantage over direct imaging for regions close 
to the star, typically over the range of 2 — 5X/D for our initial measurements presented in 
this paper. By suppressing the diffraction pattern and the quasi-static speckles that are tied 
to them, we will then be hmited by the noise characteristics of the uncorrected AO halo, 
which are only just becoming significant in our data set. Further observations will determine 
the effectiveness and sensitivity wc will reach with the APP in hour long exposures, and 
combined with a refinement in the APP models, wc will be able to detect point sources at 
much smaller separations than is possible with direct imaging and classical coronagraphs 
alone. 

The suppression of diffraction limits the quasi-static speckle noise that plagues the de- 
tection of point sources in imaging surveys for faint companions at small angular separations. 
The APP is simple, efficient and robust compared to other diffraction suppression methods, 
and can be used in conjunction with other coronagraphic designs, making it a potentially 
useful aid in the quest for direct imaging of extrasolar planets. 

This work was supported by NASA under grant NNG06GE06G and has been encouraged 
by discussions with Nick Woolf. We are grateful to Gregg Davis at II- VI, and thank Elliot 
Solheid for drawing up the filter masks and holders for the phase plate. We also thank our 
anonymous referee for the many helpful comments and suggestions which helped improve 
this paper. This work is based upon work supported by the National Aeronautics and Space 
Administration through the NASA Astrobiology Institute under Cooperative Agreement 
No. CAN-02-OSS-02 issued through the Office of Space Science. The mathematical theory 
was developed in part under grants from NASA, APRA04-0013-0056, and the NSF, AST- 
0138347. 
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Table 1. Position Angles, Separations, Differential Magnitudes and Filters of ji Her A 



Date 


P.A. 


Separation 


Am 


Filter 


2006.28070 


214.55 ±0.20 


1.4907 ±0.0048 




M 


2006.28072 


214.69 ±0.18 


1.4847 ±0.0049 


4.99 ± 0.06 


M' 



I 

to 

I 



